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ABSTRACT 


This paper examines the effect of a power law gravity field 
on barocllnlc instability. We focus on the case of inverse fifth power 
gravity, since this is the power law produced when terrestrial gravity 
is simulated in spherical geometry by electrostatic means. We have 
obtained growth rates of unstable normal modes as a function of parameters 
of the problem by solving a second order differential equation numerically. 
Results are compared with those from an earlier study where gravity was a 
constant. The conclusion is that, over the range of parameter space ex- 
plored here, there is no significant change in the character of theoretical 
regime diagrams if the vertically averaged gravity is used as parameter. 


1. Introduction 


It is generally agreed that the wavelike disturbances seen in 
rotating cylindrical geometry laboratory flow experiments are the result of 
barocllnic instability of axiaJ ly-symmetric flow. This interpretation has 
given a good theoretical basis by the study of Barcilon (196A) , in 
which the eady model of barocllnic instability was used to obtain stability 
criteria for axially-symmetrlc flows. These criteria are in reasonable 
agreement with those observed in the laboratory for rotating anmllus 
flows. 

In the earth’s atmosphere barocllnic instability is also an impor- 
tant process for maintaining departures from axially-symmetric flow. It has 
long been felt that better simulation of atmospheric flow patterns or at 
least a better understanding of how and when barocllnic Instability operates * 
on the atmosphere could be achieved if laboratory rotating fluid experiments 
could be done in spherical geometry. Such experiments have not been realize- 
able because the dielectric body force for simulated radial gravity cannot 
be made large enough to dominate the effect of ambient terrestrial gravity 
in the laboratory. The low gravity environment aboard orbiting laboratories 
such as Spacelab, to be operational in the early 1980’ s, affords an oppor- 
tunity for such an experiment. 

In going from cylindrical geometry to spherical geometry in a 
rotating fluid experiment, one important new feature is the latitudinal 
variation of the local vertical component of rotation. As is well known, 
the effect of this on the dynamics of low frequency geophysical motions 
can be taken into account by B-plane geometry. 'As one of the first steps in 
developing a model for use in design of a Spacelab experiment, Geisler and 
Fowlis (1979) extended the work of Barcilon (1964) to a g-plane channel. 

The principal result of theirstudy was to document the changes in the shape 
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and location of the baroclinlcally unstable region of parameter space 
brought about by the latitudinal dependence of the vertical component of 
rotation. 

One consequence of using a dielectric body force to simulate 
gravity is that the force field law is one of Inverse fifth power (Hart, 

1976). This must be taken into account in mathematical models of the pro- 
posed experiment and, moreover, is potentially troublesome because it does not 
simulate the inverse square of terrestrial gravity. 

This paper describes the extension of the baroclinic instability 
model of Geisler and Fowlls (1979) to Include an inverse fifth power law of 
gravity. The study shows that there is little difference between the stability 
information obtained from the two models provided gravity is replaced by its 
vertical average. This result supports conclusions earlier obtained for an 
annular geometry model by Giere and Fowlis (1979). 

2. Formulation 

Baroclinic instability in the presence of constant gravrity was 
treated in Geisler and Fowlls (1979). In that paper we obtained growth rates 
and eigenfunctions for unstable modes in both the Chamey and Eady models of 
baroclinic instability with and without Ekraan damping at the boundaries. 

The normal modes were assumed to have the functional form 

i);(x,y,z,t) = (Kz) sin ( mry \ exp [ik(x-ct)] (1) 

V W 

where is a stream function, h is the width of the channel and n is a 
positive integer. Tlie equations solved there was 
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Wlu’ro the plus sipn applies at the upper boundary and the minus sign applies 
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at the lower boundary. 

In the above equations U(z) is the basic state flow whose stability 
is being examined. As in Ccisler and Fowlls (1979), we assume that the part of 
basic state temperature field T(y,z) associated with 11 (z) decreases linearly 
with y and the basic state temperature field <T(z)> associated with the static 
stability of the fluid increases linearly with z. These paremeters enter the 
problem through the thermal wind equation. 

^ /il\ 

dz “ ■ f^ \Dy ; 
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and through the definition of N 

» go ^d<T>j (5) 

Here T is the zonally averaged temperature, <•!> is the area averaged temper-' • 
ature, a is the coefficient of thermal expansion of the fluid and g is gravity. 
In this paper we assign to g the variation. 

g “ e„ (6) 

(1 + z/a)^ 

where p is an integer, a is the inner radius of the laboratory device and is 

the value of g at z ■ 0. Integration of (4) upward from z ■ 0 (where we take 

U • 0) gives rhe basic state flow 

U(z) - _ag^a 3T F 1 
f^ 3y 

If p * 1, integration of (4) gives the logarithmic flow 



The stability of the flow given by equation (8) was examined by (liere and Fowlis 
(1979). In the present paper we examine the stability of tlie flow given by 
equation (7) with p « 5 , tliat is, an inverse fifth power gravity. 

In the case of u (z) more general tiian the linear variation with z 
used in t^Msler and Kowlis (1979), the parameter in equation (1) should 


j 1 - d+z/a)j' 
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be replaced by 


d / f ‘ cli\ 


( 9 ) 


However, as can be seen from equations (4) and (5), the factor g cancels out 
the correction to 6 then vanishes when 3T/3y and d<T>/d* are constant, as is the 
case here. 

We introduce non-dimensional quantities denoted by a prime as follows: 


x/L 

y/L 

z/d 


k' = kL 
U ’ - U /Af. 


( 10 ) 


z' ■ z/d c' ■ U /Al' 

Here L is an arbitrary horizontal length scale (taken to be 0.707a in Geisler 
and Fowlis, 1979)and dis the depth of the fluid. The quantity All is u (d) , that 
is, the difference in tlie basic state flow between z ■ 0 and z ■ d. Equations • 
(2) and (3) then become 
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Tlie parameters in these oqu at ions are a static stability S, a P-parameter B, 


an Ekman number E and a thonnal Rossby number R . Tliev are defined as 
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In the sequel, wo refer to models with E *= 0 in (12) as Ijivi.'^eid 
mi>delj>. Equations (11) and (12), we refer to as the Clun* ne^^ 
baroelinic instability. The model obtained by setting, B = 0 in (11), we 
rt'fer to as the E:u|v mot^el. 

In the invlseid Eidy imHlel it is eusttMnary to display the behaviour 
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of C| (the imaginary part of the mode phase speed) in a diagram of KS*c*^ 

L 

versus KS . Here S is given by (13) and 

k'2+nW N"* (1 

Glere and Fowlis (1979) studied the unstable modes in the inviscid Eady 
model using the log profile basic state flow given by equation (8). They 

discovered that the curve changed very little over a wide range of the 

— L • L — 

parameter d/a if S is used in the plot instead of S , where S is the value 

2 

of S obtained from using in the definition of N (see equation (6) and (13)) 
the vertically averaged value of g defined as 




In the results reported in this paper we also use S rather than S In the graphical 
presentation. 

3. Inviscid Model Results 

All results presented here and in Section A were obtained by solving 
equation (11) subject to boundary conditions (12) using the numerical technique 
briefly described in Geisler and Fowlis (1979). As noted in Section 2, the 
difference between that and the present study is that we here use gravity 
which varies with radial distance according to the function 


8 “ 8 ^ ( 17 ) 

where a is the inner radius and d the depth of the fluid. 

As a cross check we first considered the case of p » 1 in the 
inviscid Eady model, a case for which analytic solutions were obtained by 
Gierc and Fowlis (1979). Wc show in Figure (1) a plot of for runs with 
d/a = 0 (that is, constant g) , d/a = 1 , and d/a = 10 , The figure confirms 
their result that inverse first power gravity has little effect provided the 


5 


vertical average gravity is used for the plot* This agreement with their 
analytic result also indicates that the numerical routine is functioning 
correctly. 

Figure (2) shows results for the inviscid Eady model when p ■ 3 
(inverse fifth power gravity). Runs shown are those for d/a ■ 0, d/a »0.2 
and d/a * 1. We have not gone beyond d/a ■ 1 because we anticipate that d/a 
will have to be <1 in any reasonable geophysical experiment. Figure (3) shows 
results when we go to the inviscid Chamey model. Both Figure (2) and Figure 
(3) support the conclusion that, for the range of d/a of geophysical interest 
at least, such a plot is little affected by inverse fifth power gravity if the 
vertical average gravity is used. 

4. Regime Diagrams 

Superimposed contour diagrams of c^ for several values of zonal wave 

2 '' 

number k in the parameter space of S versus R^ /S'"E constitute a theoretical 

r egime diagram . The construction of these diagrams for the Eady model and the 

Chamey model was the subject of the paper by Geisler and Fowlis (1979). 

Figure (4) is taken from that paper and shows a theoretical regime diagram for 

2 2 2 

the Eady model when d « h • L - 0.707a. (The parameter e /6 « R^ /S E). The 
curve labelled (1) is the stability bounaary for zonal wave number (1) disturb- 
ances, outside (to the left) of which c^ <0. The envelope of the curves shown 
in this figure obviously separates the region of parameter space where the flow 
is unstable from where the flow is stable and hence axially-symmetric. Geophysical 
experiments must be such that the unstable regime can exist in the apparatus, 
hence the utility of theoretical regime diagrams in experimental design studies. 

To illustrate the results of our study of the effect of Inverse fifth 
power gravity on theoretical regime diagrams, we have selected zonal wave 
number (3). Figure (5) shows contours :>f KS for this wave number in the 
Eady model when g is constant. The second case, where d/a « 0.5, is so close 

to this that it does not bear showing. Figure (6) shows the extreme (for a 
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geophysical experiment) case for d/a ■ 1. Comparison with Figure 5 shows 
that the difference is rather smoll. We have run cases for other zonal 
wave numbers with the same result. The conclusion is that for d/a ^1, 
the inverse fifth power gravity has only a small effect on the shape and 
location of the unstable regime provided vertically averaged g is used in 
drawing the diagram. 

— *5 

Figure (7) shows contours of KS^c^^' for zonal wave number (3) in 
the Chamey model when g is constant. Here we adopt the value B ■ 2.35 used 
as a standard case in Geisler and Fowlis (1979). The case for d/a - 0.2 is 
again not much different and is not shown here. Figure (8) shows the case for 
d/a * 1. The change is somewhat greater than in the corresponding Eady model 
runs, the most notable change occurring in the region S _^0.l. However, there 
is very little change in the leftward penetration of the nose located at about 
S 0.2. We have obtained similar results for other zonal wave numbers. We 
conclude that for d/a the Inverse fifth power gravity does not have significant 
affect on regime diagrams in either the Chamey model or the Eady model provided 
S rather than S is used as the parameter. 
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Fisure Letends 


Figure 1. Inaginary part of phase speed in the Inviscid Eady node! for the 
case of Inverse first power gravity. Curves are labelled by the 
value of d/a. The dashed curve is the case of constant gravity. 

Figure 2. Inaginary part of phase speed in the inviscid Eady model for the 
caae of inverse fifth power gravity. Curves are labelled by the 
value of d/a. The dashed curve is the case of constant gravity. 

Figure 3. Same as for Figure 2, but in the case of the inviscid Chamey 
model. 

Figure 4. Theoretical regime diagram for the Eady model with d ■ h •• acos6^ 
(from Gelsler and Fowlis^ 1979). The curves are atability 
botsidaries labelled by zonal wave number. 

Figure 5. Contours of imaginary part of phase speed KS c^* in the Eady model 
with gravity constant. . 

Figure 6. Contours of imaginary part of phase speed KE’c^' in the Eady model 
with Inverse fifth power gravity and d/a ■ 1. 

Figure 7. Contoura of Imaginary part of phase speed KS^c. * in the Chamey 
model with gravity constant. 

Figure 8. Contours of imaginary part of phase speed KS i" the Chamey 
model with inverse fifth power gravity and d/a *• 1. 
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